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non-Gaussian  noise  is  discussed.  The  optimal  detector  non¬ 
linearity  is  approximated  adaptively  in  the  noise  pdf  tail  region, 
and  a  polynomial  is  used  to  approximate  the  non-linearity  near  the 
mean.  Examples  for  several  different  noise  environments  are  pre¬ 
sented,  showing  in  these  cases  that  the  adaptive  detector  is  able 
to  achieve  a  high  level  of  performance 


ADAPTIVE  DETECTION  OF  A  KNOWN  SIGNAL  IN  NON-GAUSSLAN  NOISE 


STEVEN  V.  CZARNECM  and  JOHN  B.  THOMAS 
Department  of  Electrical  Engineering  and  Computer  Science 
Princeton  University 
Princeton,  New  Jersey  08540 


ABSTRACT 


•  -  .4  (* 

;  «  «  »  *-a 


Hie  design  of  a  locally  optimal  detector  for  a  known  signal  in  non-Gaussian 
noise  is  discussed.  Hie  optimal  detector  nonlinearity  is  approximated  adaptively 
in  the  noise  pdf  tail  region,  and  a  polynomial  is  used  to  approximate  the  non¬ 
linearity  near  the  mean.  Examples  for  several  different  noise  environments  are 
presented,  showing  in  these  eases  that  the  adaptive  detector  is  able  to  achieve  a 
high  level  of  performance 
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1.  INTRODUCTION 

A  binary  hypothesis  test  may  be  used  to  model  the  problem  of  detecting  a 

known  signal  in  the  presence  of  noise.  Let  ds  =  |sj . be  the  known  signal 

sequence  with  amplitude  parameter  d  >  0,  and  let  n  =  jnj . ny\  be  an  iid 

noise  sequence  independent  of  the  signal.  The  detector  observes  x,  a  data 
sequence  [xy  ...  ,xy\,  and  decides  between: 

H0.  x  =  n 
Hy  x  =  n+ds 


In  the  framework  of  Neyman-Pearson  hypothesis  testing,  x  and  the  multivariate 
noise  density  /  are  used  to  calculate  a  likelihood  ratio  ANP.  This  test  statistic 
and  a  fixed  threshold  TNP  are  compared  to  arrive  at  a  decision:  is  chosen 

when  Ajvp  >  TNP,  and  HQ  is  chosen  when  AyP  *  Typ.  More  precisely, 


_  /ifrlffi)  _  /x(x-d s) 

'  /kW*o>  '/*(*)  *  NP 


(I) 


Since  the  noise  is  iid,  and  the  logarithm  function  is  monotonic,  an  equivalent  test 
is 


\ NP  ~  111 A HP  ~  S 9NP,i\xi)  -  ^  tffP  =  hlTyp 

»■!  <*1  ^ 


(2) 


where  /n  is  the  univariate  density  of  n*.  When  the  signal  is  constant,  s,  =  s  for 
*  «  L  —  ,Af,  and  the  sequence  {gyp*  {  may  be  replaced  with  a  zero  memory 
nonlinearity,  gyP. 

In  cases  where  the  aignal-to-noise  ratio  is  very  small,  the  test  statistic  may 
be  calculated  via  the  locally  optimal  detector.  The  test  becomes 


*10  -  £  VLo(*i)  =  "  ^  *  *LO 


i  •  1 


t  ■  r 


(3) 


U  o 


Implicit  in  both  detection  methods  is  a  requirement  that  the  noise  pdf  must  be 
known  exactly.  In  general,  the  noise  statistics  are  not  known  with  precision  and 
the  design  of  the  LO  or  NP  detector  is  not  straightforward.  Alternative  detection 
strategies  are  available,  and  among  these  are  (l)  detectors  which  are  robust  with 
respect  to  deviations  from  a  nominal  noise  environment  [2,8];  (2)  nonparametric 
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delectors  which  use  only  very  general  information  about  the  underlying  noise 
distribution  [3];  and  (3)  fixed  or  adaptive  suboptimal  detectors  which  have  nearly 
optimal  performance  [4.7,13]-  An  additional  consideration  is  that  often  the  noise 
environment  is  nonstationary  and  an  adaptive  structure  is  necessary. 

2.  DETECTION  IN  NON-GAUSSIAN  NOISE 

The  classical  assumption  in  the  design  and  analysis  of  detection  systems  is 
that  the  noise  is  Gaussian.  This  assumption  has  attractive  features  in  that  the 
Gaussian  model  is  mathematically  tractable,  the  optimal  detector  structure  is 
linear,  and  strong  justification  for  the  model  is  available  in  the  form  of  the  Cen¬ 
tral  limit  Theorem.  Measurements  of  typical  noise  environments  have  led  to  the 
conclusion  that  the  true  noise  distribution  often  is  described  better  by  a  heavier 
tailed  pdf  [5,6,12,13].  This  type  of  noise  can  be  ascribed  to  a  nominal  Gaussian 
environment  with  a  heavy  tailed  impulsive  noise  contaminant,  and  to  the  fact 
that  a  real  noise  comprises  a  finite  sum  of  random  events;  convergence  to  the 
Gaussian  pdf  is  not  complete  for  a  finite  sum,  but  the  sum  pdf  is  most  nearly 
Gaussian  near  the  mean,  with  the  tails  converging  to  the  Gaussian  pdf  only  in  the 
limiting  case.  For  example,  previous  work  in  detection  has  taken  note  of  non- 
Gaussian  noise  environments,  where  the  noise  density  is  approximately  Gaussian 
near  the  mean,  but  heavier  tailed  away  from  the  mean  [1,6,10], 

The  NP  and  LO  detector  nonlinearities  related  to  near-Gaussian  heavy  tailed 
densities  are  typically  composed  of  a  linear  region  surrounded  by  tails  which 
compress,  limit,  or  even  blank  large  data  observations.  Work  in  robust  estima¬ 
tion  of  the  mean  has  similarly  suggested  that,  in  heavy  tailed  noise,  a  robust 
estimator  should  reduce  the  influence  of  very  large  data  observations  while  leav¬ 
ing  observations  near  the  mean  relatively  unchanged  [9]. 

3.  A  SUBOPTIMAL  DETECTOR  APPROXIMATION 

This  paper  presents  an  approach  to  the  design  of  a  noise-adaptive  subop¬ 
timal  detector  with  these  ideas  in  mind.  Attention  is  focused  on  the  case  of  LO 
detection  of  a  constant  binary  signal  in  discrete  time,  with  iid  noise  assumed. 
Further,  the  noise  pdf  is  restricted  to  be  unimodal,  symmetric  about  its  mean 
placed  at  the  origin,  and  to  have  nonzero  support  over  the  entire  real  line.  We 
are  interested  in  approximating  ZNL's  loosely  specified  by  the  following  charac¬ 
teristics: 

a)  continuous,  with  continuous  low-order  derivatives 

b)  approximately  linear  at  the  origin 

c)  odd  symmetric  about  the  origin 

d)  strictly  positive  to  the  right  of  the  origin 

e)  monotone  in  the  tail  regions 

f)  possess  only  one  local  extremum  on  either  aide  of  the  origin 

Pdf's  which  are  non-Gaussian  in  the  sense  previously  described  can  generate  LO 
nonlinearities  with  these  characteristics,  which  are  typical  of  robust  and  subop¬ 
timal  detector  nonlinearities,  and  also  of  many  of  the  influence  curves  given  for 
robust  estimation  procedures  [  1 1]. 

The  tail  behavior  of  these  LO  nonlinearities  range  from  linear  for  a  noise  pdf 
with  Gaussian  tails,  to  a  limiter  for  exponentially  decreasing  pdf  tails,  to  a 
blanker  for  algebraically  decreasing  pdf  tails,  bi  general,  the  heavier  tailed  the 
noise  density  is  relative  to  the  Gaussian  pdf,  the  more  severely  curtailed  is  the 
effect  of  large  data  observations. 

The  objective  of  a  noise  adaptive  nonlinearity,  then,  should  be  to  relate  the 
ZNL  tail  behavior  to  the  actually  observed  noise  pdf  tail  behavior.  We  propose 
the  following  method:  It  has  been  reported  [12]  that  the  generalized  Gaussian 
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density 

/e(»)  =  2ru/cTexp(~|T?g|c^ 

in  certain  instances  can  describe  the  pdl  tails  o!  some  physical  noise  sources. 
For  a  noise  variance  of  a2,  the  parameter  rj  is  defined  by 

V 

The  corresponding  LO  nonlinearities,  shown  in  Figure  1,  can  be  written  as 

9lo(*)  =  crf  I*  !e-,sgn(x)  (5) 

with  the  c  conveniently  parameterizing  ZNL  tail  behavior.  Therefore,  we  model 
the  observed  noise  pdf  tails  via  the  generalized  Gaussian  family.  The  suboptimal 
LO  nonlinearity  thus  will  have  power  law  tails  described  by 

§lo(*)  =  Srj'lxl^sgnfx)  for|x|>x0  (6) 

It  is  necessary  to  find  a  value  c  such  that  /e  is  a  good  approximant  to  the  tail 
behavior  of  /n ,  the  true  but  unknown  underlying  noise  density.  A  simple  way  to 
do  this  is  to  equate  the  tail  probability  of  ft  with  the  observed  tail  mass 

Pt  ~  Ar.-)dr*  I)  (7) 

ni  ■  i 

Here,  I  is  the  indicator  function  and  are  the  noise  observations  presumed 
available  from  a  noise  reference  channel.  The  exponent  c  can  be  estimated  then 
by  finding  8  such  that 

aa 

2  ffe(x)dx  =  PT  (8) 

T 

For  the  purpose  of  calculating  rj,  the  observed  noise  variance  a  at  2  is  used.  The 
estimate  c  is  defined  by  the  implicit  integral  in  (B),  so  it  is  desirable  to  derive  a 
simpler  explicit  relationship 

8  =  hf (Pf)  (9) 

One  obvious  method  is  to  first  calculate  Pf  as  a  function  of  c ,  and  then  interpo¬ 
late  this  tabulation  to  approximate  (9).  This  tabulated  version  of  Ay  is  shown  in 
Figure  8. 

With  a2  fixed  and  c  small,  the  value  of  r?,  a  scale  factor,  becomes  large. 
Therefore,  even  though  /c  (x)  approaches  zero  asymptotically  at  a  much  slower 
rate  than  the  Gaussian  pdf  as  |x  j  gets  large,  the  total  probability  mass  in  the 
tails  is  quite  small.  As  a  result,  hT  is  multiple  valued,  the  density  is  extremely 
peaked,  and  the  LO  nonlinearity  has  a  discontinuity  at  the  origin.  Clearly,  for 
c  <  1,  the  requirement  of  near-linearity  at  the  origin  is  not  met  by  (S).  Since  the 
objective  in  using  the  generalized  Gaussian  pdf  is  to  relate  the  tall  heaviness  of 
an  observed  noise  to  a  parameter  governing  the  shape  of  the  ZNL  tail,  we  replace 
the  anomalous  behavior  of  the  true  function  A  with  a  simple  linear  relation 

&r(£r)  *  oPr+b  (10) 

where  a  and  6  are  chosen  to  approximate  (B)  for  a  particular  value  of  T.  This 
approximation  is  plotted  as  the  broken  line  on  Figure  2.  The  values  for  a  and  6 
are  chosen  so  that  when  Pf  corresponds  to  Gaussian  or  exponential  noise  tails, 
(10)  gives  8=8  and  8sl,  respectively,  fcote  that  the  linear  relation  allows  the 
value  of  8  to  be  negative  for  large  tail  probabilities. 


The  tail  measurement  threshold  T  must  be  chosen  prior  to  estimating 
parameters  a  and  6 .  One  way  to  pick  T  is  to  choose  a  value  which  minimizes 

Ep  jrarej  =  a%  — UD 

where  N  is  the  number  of  noise  observations.  For  c  uniformly  distributed  on  the 
interval  [1,2),  and  Pj  the  tail  probability  of  /e (*)  lor  |z|  >  T,  the  constant 
T  -3a  approximately  minimizes  (11). 

The  LO  nonlinearities  of  the  generalized  Gaussian  family  have  desirable  tail 
behavior,  but  for  small  values  of  £  the  behavior  does  not  meet  the  constraint  of 
linearity  near  the  origin.  To  eliminate  this  behavior,  the  ZNL  needs  modification 
in  the  region  near  the  origin.  A  way  in  which  to  do  this  is  to  replace  §io(x),  for  z 
near  zero,  with  a  function  that  fits  the  given  desired  characteristics  (a-f).  A  suit¬ 
able  family  of  functions  are  polynomials  p(x)  with  the  following  characteristics: 

p(x)  =  a3z8+a2z8+a1z+o0  for  0£z£z0 

P(0)  =  0 

J>(l±*ol)s£n(±*o)  =  5io(±*o) 

P'(l±*o!)  =  0'.w<±*o) 
p*'(l±z0|)sgn(±z0)  =  $"io(±z0) 

The  choice  of  tail  behavior  via  c  and  the  point  zo  completely  specify  p(x).  The 
method  for  choosing  c  has  already  been  specified,  so  the  choice  of  z0  is  the 
remaining  free  parameter. 

A  method  equivalent  to  choosing  the  proper  z0  is  to  choose  an  arbitrary  z0 
and  scale  the  input  to  the  ZNL  with  a  factor  v.  It  is  reasonable  to  choose  v  to 
maximize  the  efficacy  of  the  ZNL.  For  an  arbitrary  nonlinearity  q,  efficacy  as  a 
function  of  v  can  be  defined  as 

i  („)  _ _ ^E/[j'(w)] _ 

4,1  ’  1  ’ 

where  E/  refers  to  the  expectation  with  respect  to  the  density  function  /(z), 
and  q'  refers  to  the  first  derivative  of  q . 

At  this  point,  specification  of  the  suboptimal  nonlinearity  g  is  complete,  and 
can  be  written  as 

v  p(|uz|)sgn(vz)  **  l«l**o 

=  l  B|»|’->«gn(i«)  If  Iwtl  >*o  (13) 

Figures  3a-3c  gives  some  examples  of  the  types  of  nonlinearities  available  using 
this  approximation  method. 

4.  EXAMPLES 

Analytic 

We  will  now  present  examples  of  the  use  of  y  in  approximating  some  known 
optimal  LO  nonlinearities. 

The  first  comparison  is  between  the  approximate  and  exact  versions  of  LO 
nonlinearities  for  the  generalized  Gaussian  family.  The  exponent  S  is  given  by 
(10).  after  using  the  exact  value  c  in  (8)  to  obtain  Pf.  Since  this  is  an  analytical 
example,  and  the  true  noise  density  is  known,  numerical  methods  can  be  used  to 
obtain  v*.  the  value  of  u  which  maximizes  £v(».  The  performance  of  the  subop¬ 
timal  ZNL  relative  to  the  LO  nonlinearity  can  be  measured  by  asymptotic  relative 
efficiency,  which  is  given  easily  in  terms  of  efficacy  as 


Figure  4  compares  the  performance  of  g,  the  LO  detector  and  a  linear  detector 
(Id),  In  terms  of  AHEJjio  and  AREgwM.  The  suboptimal  nonlinearity  performs 

quite  well  for  the  range  U  c  s  2,  but  for  c  <  1,  performance  deteriorates.  This 
is  easily  explainable,  since  for  small  c ,  the  LO  nonlinearity  output  approaches  ±°° 
for  inputs  near  zero,  while  the  approximation  method  requires  g  to  pass  through 
the  origin 

Another  family  of  heavy  tailed  pdf’s  is  the  Johnson  Su  family.  If  A'  is  distri¬ 
buted  as  an  N(0,<r)  random  variable,  and  we  define  a  new  random  variable 


then  the  density  of  Y  has  unity  variance,  and  belongs  to  the  Johnson  Su  family, 
given  by 

Uv)  =  x&rflx]  <i6> 

with 


The  parameter  6  controls  the  tail  heaviness.  As  <5  -*  <*>,  the  pdf  tails  become  pro¬ 
gressively  lighter,  and  approach  Gaussian  tails  in  the  limit. 

Since  /  a  is  given  and  known.  Pf  can  be  calculated  from  (0),  and  (10)  gives  £. 
Again,  numerical  methods  can  find  the  v*  which  maximizes  (12).  Some 
representative  LO  nonlinearities  and  suboptimal  approximations  are  given  for 
various  values  of  6  in  Figure  5,  and  Figure  6  presents  the  performance  com¬ 
parison  of  g,  gio,  and  Id.  For  this  set  of  densities,  the  approximation  method 
works  quite  well.  Over  the  range  .0  £  6  ss  »,  the  minimum  of  AREj^  is  .909, 

(occurring  for  6  ~  .0).  This  means  that  only  a  small  performance  penalty  would 
be  incurred  if  g  were  to  replace  the  LO  detector.  As  a  final  comment,  it  should 
be  observed  that  unlike  the  generalized  Gaussian  family,  the  Johnson  Su  family 
fulfills  the  characteristics  of  a  near  Gaussian  pdf  given  earlier. 

A  third  family  of  heavy  tailed  densities  was  also  investigated.  This  third  fam¬ 
ily  is  the  e -contaminated  Gaussian  mixture,  and  in  this  case,  the  Gaussian- 
Gaussian  mixture  The  pdf  of  this  noise  can  be  written  as 

/«(*)  «  (l-e)/0(*)  +  e/i(*)  (18) 

where  /0  represents  the  pdf  of  an  N(0,1)  distributed  random  variable,  and  /] 
represents  the  pdf  of  an  N(0,o*)  random  variable,  with  o*  large.  The  parameter 
e  controls  the  degree  to  which  /  j  contaminates  the  nominal  density  /o.  and  is 
typically  taken  to  be  small.  The  performance  of  the  suboptimal  ZNL  was  calcu¬ 
lated  for  some  specific  cases,  and  the  results  appear  in  Table  1.  Figure  7  shows 
the  comparison  in  the  form  of  two  LO  nonlinearities  and  the  corresponding 
approximations  The  performance  of  the  approximations  in  this  case  is  not  as 
good  as  the  approximations  under  the  Johnson  5^  family.  Once  again  not  all  the 
assumptions  are  met  by  the  c-contaminated  Caussian-Gaussian  mixture  on  the 
form  of  the  LO  nonlinearity,  as  there  is  more  than  one  local  extremum  on  each 
side  of  the  origin. 


-6 


Table  1 

Performance  of  gu>  and  g  in  c-contaminated  Gaussian-Gaussian  noise 


e 

*0 

c 

u* 

AREJi0,u 

20 

.05 

3 

.734 

.767 

.952 

1.72 

20 

.20 

3 

-.196 

.821 

.960 

3.10 

5 

.05 

3 

1.54 

.957 

.977 

1.08 

5 

.20 

3 

.996 

.957 

.983 

1.24 

The  results  so  far  show  that,  in  several  analytical  examples,  it  is  possible  to 
achieve  reasonably  good  performance  from  this  simple  approximation. 
Simulation 

To  see  how  well  this  system  might  work  in  practice,  some  actual  physical 
noise  was  used  to  drive  the  system.  The  noise  was  collected  underneath  the  Arc¬ 
tic  ice  pack,  and  details  may  be  found  in  [13],  The  noise  data  is  highly  nonsta¬ 
tionary;  a  background  Gaussian  noise  is  abruptly  interrupted  with  segments  of  a 
very  heavy  tailed  noise  generated  by  cracking  of  the  Arctic  ice  pack.  The  kur- 
tosis  for  each  consecutive  block  of  1024  noise  samples  was  calculated,  and  blocks 
with  a  kurtosis  exceeding  4  were  selected.  To  get  a  more  nearly  stationary  noise 
for  driving  the  system,  the  data  in  each  block  was  normalized  to  unit  variance, 
randomly  permuted,  and  the  blocks  were  then  concatenated,  thus  simulating  the 
output  of  a  physical  noise  source.  This  noise  was  used  as  the  input  for  the  adap¬ 
tive  system.  The  exponent  c  was  estimated  from  the  running  estimate  of  Pf, 
which  converges  to  the  true  tail  probability  as  N -*<*>.  The  scale  parameter  v  can¬ 
not  be  found  via  numerical  methods,  since  in  this  case  the  simulated  system  has 
no  knowledge  of  the  true  density  generating  the  noise  observations.  However, 
the  fGefer-Wolfowitz  stochastic  approximation  method  can  be  used  to  find  the 
value  of  v*  which  maximizes  The  convergence  rote  towards  v*  is  fixed  by  the 
particulars  of  the  SA  algorithm,  and  no  formal  attempts  were  made  to  optimize 
its  performance. 

Figure  8  shows  the  running  estimate  of  c  and  v  as  a  function  of  sample 
number.  Figure  9  shows  the  estimated  efficacy  of  g  for  each  1024  sample  block, 
and  the  cumulative  estimate  of  efficacy  starting  from  the  first  noise  observation. 
At  the  end  of  the  simulation,  it  was  assumed  that  the  approximated  ZNL  was  as 
near  optimal  as  possible.  The  final  ZNL  is  shown  in  Figure  10.  The  efficacy  of  this 
suboptimal  detector  was  calculated  using  the  empirical  distribution  of  the  entire 
noise  observation  set,  and  it  was  found  to  be  1.39.  The  true  distribution  of  the 
noise  is  not  known,  and  therefore  it  is  not  possible  to  calculate  the  efficacy  of  the 
true  L0  nonlinearity.  However,  it  is  possible  to  conclude  that  the  adaptive  sys¬ 
tem  was  able  to  use  the  noise  observations  and  adapt  the  detector  structure  in  a 
constructive  way.  The  ARE  of  the  approximate  ZNL  relative  to  the  linear  detec¬ 
tor  is  given  simply  by  the  efficacy  in  this  case,  and  therefore  the  adaptive  ZNL 
shows  an  improved  performance  over  the  linear  detector. 


5.  CONCLUSION 

The  conclusion  to  be  drawn  from  this  study  is  that  it  is  possible  to  imple¬ 
ment  an  adaptive  detector  nonlinearity  using  fairly  simple  techniques.  The  esti¬ 
mate  of  tail  behavior  is  quite  simple:  merely  a  measurement  of  the  relative 
number  of  samples  exceeding  a  specified  threshold.  This  apparently  gives 
enough  information  about  tail  behavior  of  the  true  noise  density,  so  that  even  a 
crude  approximation  to  the  true  ZNL  in  the  tail  regions  results  in  fairly  good  per¬ 
formance.  A  more  sophisticated  estimate  of  C  might  improve  the  performance 
g  It  would  be  interesting  to  discover  how  much  additional  complexity  any 
resulting  performance  gain  could  justify. 
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Some  natural  amplifications  to  this  approximation  method  suggest  them* 
selves  for  investigation.  One  is  to  replace  the  central  region  polynomial  in  the 
ZNL  simply  by  a  linear  connection  between  the  tail  approximants  at  ±x0.  A 
further  simplification  would  be  to  replace  the  power-law  tails  with  piecewise 
linear  approximations,  possibly  even  a  single  straight  line.  Some  work  done  by 
others  [7J  in  approximation  of  LO  nonlinearities  suggest  that  even  very  simple 
approximants  of  the  optimal  nonlinearity  have  the  potential  to  achieve  perfor¬ 
mance  which  is  acceptably  near  the  optimal. 

ACKNOWLEDGEMENT 

This  research  is  supported  by  the  Office  of  Naval  Research  under  Contract 
NOOO 14-8 1-K-0 146  and  by  the  National  Science  Foundation  under  Grant  ECS- 79- 
10915. 

REFERENCES 

[1]  J.  H.  Miller  and  J.  B.  Thomas,  "Detectors  for  Discrete-Time  Signals  in  Non- 
Gaussian  Noise,"  IEEE  Trans.  Inform.  Theory,  vol.  IT-18,  no.  2,  pp.  241-250, 
March  1972. 

[2]  P.  J.  Huber,  "A  Robust  Version  of  the  Probability  Ratio  Test,"  Ann.  Math.  Sta¬ 
tist.,  vol.  36,  pp.  1753-1758,  Dec.  1965. 

[3]  J.  B.  Thomas,  "Nonparametric  Detection,"  Proc.  IEEE,  vol.58,  no.  5,  pp.623- 
631,  May  1970. 

[4]  J.  W.  Modestino,  "Adaptive  Detection  of  Signals  in  Impulsive  Noise  Environ¬ 
ments,"  IEEE  Trans.  Commun.,  vol.  COM-25,  no.  9,  pp.  1022-1026,  Sept.  1977. 

[5]  S.  L.  Bernstein,  et  at.  "Long-Range  Communications  at  Extremely  Low  Fre¬ 
quencies,"  Pro c.  IEEE,  vol.  62.  no.  3,  pp.  292-312,  March  1974. 

[6]  A.  D.  Spaulding  and  D.  Middleton,  "Optimum  Reception  in  an  Impulsive 
Interference  Environment  -  Part  I:  Coherent  Detection,"  IEEE  Trans.  Com¬ 
mon.,  vol  COM-25,  no.  9,  pp.  924-934,  Sept.  1977. 

[7]  J.  H.  Miller  and  J.  B.  Thomas,  "Robust  Detectors  for  Signals  in  Non-Gaussian 
Noise. ”  IEEE  Trans.  Oommun.,  vol.  COM-25,  no.  7,  pp.  6B6-690,  July  1977. 

[8]  A.  H.  El-Sawy  and  V.  D.  Vandelinde,  "Robust  Detection  of  Known  Signals," 
IEEE  Trans.  Inform.  Theory,  vol.  IT-23,  no.  6,  pp.  722-727,  Nov.  1977. 

[9]  F.  R.  Hampel.  "The  Influence  Curve  and  its  Role  in  Robust  Estimation,"  Jour. 
Amer.  Stat.  Assoc.,  vol.  69,  no.  346,  pp.  383-393,  June  1974. 

[10]  R.  F.  Ingram  and  R.  Houle,  "Performance  of  the  Optimum  and  Several  Subop¬ 
timum  Receivers  for  Threshold  Detection  of  Known  Signals  in  Additive.  White, 
Non-Gaussian  Noise,"  Technical  Report  6339,  Naval  Underwater  Systems 
Center,  New  London,  CT,  Nov.  24,  1980. 

[11] D.  F.  Andrews,  st  a l,  Robust  Estimates  of  Location,  Princeton  University 
Press.  Princeton,  NJ,  1972. 

[12]  A.  D.  Watt  and  E.  L.  Maxwell,  "Measured  Statistical  Characteristics  of  VLF 
Atmospheric  Noise,"  Proc.  Inst.  Radio  Engineers,  vol.  45,  no.  1,  pp  55-62, 
Jan.  1957. 

[13] R.  F.  Dwyer,  "Arctic  Ambient  Noise  Statistical  Measurement  Results  and 
Their  Implications  to  Sonar  Performance  Improvements,"  Technical  Report 
6739,  Naval  Underwater  Systems  Center,  New  London,  CT,  May  5,  1962. 


Figure  l.  Locally  optima]  nonlinearities  for  the  generalised  Gaus¬ 
sian  family. 


Figure  2_  The  exact  relation  e  *  hT(PT)  (solid  line),  and  it*  approx¬ 
imation  e  *  ^j(Pt)  (broken  line)  for  various  T 


figure  3a.  The  ZNL  }  for  xB  «  3,  u  ■  l.  various  e . 
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Figure  3c  The  ZNL$  for  Xq  *  3.  various  v.e  *  1 


Figure  4.  Performance  of  g  and  g^j  versus  c  for  the  generalized 
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Figure  8.  $u>  (solid  line)  and  f  (broken  line)  for  Johnson  SL  family, 
on  different  vertical  scales  for  comparison  For  f*l.  c*  .782. 
v*  «  3  26,  s0  a  3  For  d  *  2,  f  *  l  48.  vm  ■  1.88,  s0  *  3. 


Figure  8.  Performance  rf  |  and  versus  5  for  the  Johnson  S, 
family.  ARE|f|f.  right  sea*.  all  others  on  left  seals 


figure  7.  fu>  (■olid  line)  and  §  (broken  line)  for  Gaussian- Gaussian 
mixture,  on  different  vertical  acalee  for  comparison  (a) 
c  *  .06  of  ■  S .  (b)  t  *  .20  ef  =  20 


figure  0.  ZNL  parameters  versus  observation  number  for  heavy 
tailed  Arctic  ice  noise 


figure  t.  Estimated  performance  of  §  for  the  heavy  tailed  Arctic 
ice  noise 


